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The Sivers function is extracted from HERMES data on single spin asymmetries in semi-inclusive 
deeply inelastic scattering. Our analysis use a simple Gaussian model for the distribution of trans- 
verse parton momenta, together with the flavor dependence given by the leading 1 /N c approximation 
and a neglect of the Sivers antiquark distribution. We find that within the errors of the data these 
approximations are sufficient. 
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I. INTRODUCTION 



Single spin asymmetries (SSA) in hard reactions have a long history dating back to the 1970s when significant 
polarizations of A-hyperons in collisions of unpolarizcd hadrons were observed [1], and to the early 1990s when large 
asymmetries in p*p — > ttX or p^p — ► 7rA were found at FNAL [2]. No fully consistent and satisfactory unifying 
, approach to the theoretical description of these observations has been found so far — see the reviews [3] . 

Interestingly, the most recently observed SSA phenomena, namely those in semi-inclusive deeply inelastic scattering 
(SIDIS) [4-11], seem better under control. This is in particular the case for the transverse target SSA observed at 
HERMES and COMPASS [8-11]. On the basis of a generalized factorization approach in which transverse parton 
momenta are taken into account [12-14] these "leading twist" asymmetries can be explained [15] in terms of the Sivers 
[16-19] or Collins effect [20]. The former describes, loosely speaking, the distribution of unpolarized partons in a 
^\ " transversely polarized proton, the latter describes the fragmentation of transversely polarized partons into unpolarized 
hadrons. In the transverse target SSA these effects can be distinguished by the different azimuthal angle distribution 
of the produced hadrons: Sivers effect cx sin(0 — <frs), while Collins effect oc sin(</> + <f>s), where <j) and tfis denote 
respectively the azimuthal angles of the produced hadron and the target polarization vector with respect to the axis 
■ defined by the hard virtual photon [15]. Both effects have been subject to intensive phenomenological studies in 
p { hadron-hadron-collisions [21-25] and in SIDIS [26-35]. For the longitudinal target SSA in SIDIS, which were observed 
JL , first [4-6] but are dominated by subleading- twist effects [36] , the situation is less clear and their description (presuming 
' the factorization theorems [12-14] can be generalized to twist-3) is more involved [37, 38]. 

In this work we will concentrate on the Sivers effect, which is quantified by the "Sivers function" fyf{x,p^) 
(in the notation recommended in [39]). It is referred to as "naively" or "artificially time-reversal-odd" for it arises 
from a correlation between the nucleon spin St and the intrinsic transverse parton momentum px, both transverse 
with respect to the nucleon momentum Pjv in the infinite momentum frame, with the effect being proportional to 
(S T x p T ) • P N - 

By adequately weighting the events entering the spin asymmetry with sin(0 — (f>s) one can project out from the data 
the Sivers SSA. Including into the weight in addition to that a power of the transverse momentum of the produced 
hadron Ph± = \Ph±\ yields an SSA which is described (with the neglect of soft factors [12-14]) model-independently 

in terms of the "transverse moment" f^? ( x ) °f the Sivers function [15]. Preliminary HERMES data analyzed 
in this way were presented [8] and subject to first studies [32]. However, the currently available final HERMES 
and COMPASS data [9, 10] were analyzed without a transverse momentum weight, and can only be interpreted 
by resorting to some model for the distribution of the transverse parton momenta in the "unintegrated" [40] Sivers 
distribution and unpolarized fragmentation function. Different models have been explored in literature [33-35]. Here 
we approximate the distribution of transverse parton momenta in the Sivers function to be Gaussian. 

We pay particular attention to the demonstration of the phenomenological consistency of the approach, and fix or 
constrain the free parameters in the Gaussian ansatz consistently by the SIDIS HERMES data. Although hereby the 
Gaussian width of the Sivers function remains poorly constrained, this does not prevent a meaningful extraction of 
the transverse moment of the Sivers distribution function from the data [9]. This demonstrates that — within the 
accuracy of the present data — the Gaussian ansatz is robust and reliable. A comparison to extractions of the Sivers 
function, where no [32] or other models [33-35] were assumed, helps to estimate the effects of model-dependence. We 
find them smaller than the statistical accuracy of the present data. 

In order to reduce the number of free parameters in the ansatz for the Sivers function, we impose as an additional 
theoretical constraint the predictions for the Sivers function from the QCD limit of a large number of colours N c [41], 
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which state that f^jf = —f^ up to l/iV c -corrections. Since a fit constrained in this way works and describes the 
HERMES data [9], our study, as a byproduct, not only tests the large- N c results [41] but also naturally explains the 
smallness of the Sivers effect from a deuteron target observed by COMPASS [10]. Besides choosing different models for 
transverse parton momenta and/or ways to fix the respective parameters, the explicit use of the large- iV c constraints 
is the main difference of our approach compared to the works [33-35] . 

For the fit we use only the HERMES data [9] on the x-dependence of the Sivers SSA. Thus, the HERMES data 
on the z-dependence serve as cross and consistency checks for the fit and the Gaussian ansatz. We also explicitly 
address the questions what could be the effects of l/jV c -corrections and Sivers antiquarks (which we neglect in our 
fit). Finally, we suggest how the Gaussian ansatz could be further tested by means of SIDIS data. Such tests are of 
importance for they allow to understand the range of applicability and the limitations of this ansatz. 



II. SIVERS EFFECT IN SIDIS 



Consider the process — > I'hX, where denotes the transverse (with respect to the beam) target polarization. 
Let P, I (?') and Ph denote respectively the momentum of the target proton, incoming (outgoing) lepton and produced 
hadron. The relevant kinematic variables are q = I — I' with Q 2 = —q 2 , x — Q 2 /(2P ■ q) and z = P ■ Ph/P ■ q- The 
Sivers SSA as presented in Ref. [9] is defined as sum over events i as follows: 



A 



sin(>-</> s ) 
UT 



<f>s,i) [A^ T (^; <t>s,i) - N^Ma^ + tt)] 



IEJATT^^+A^^ + tt)] 



(1) 



where N^^{4n;<j)s,i) are the event counts for the respective target polarization (corrected for depolarization effects) 
— see Fig. 1 for the definition of kinematics. It is understood that if the SSA is considered as function of one kinematic 
variable, then an appropriate averaging over the other variables is implied. 

In order to describe the Sivers SSA as defined in (1) we will make two major simplifications. The first consists in 
neglecting the soft factors [12-14], and their associated energy dependence. This step simplifies the description of the 
process considerably, though it is difficult to quantify the uncertainty we introduce in this way. Then, to leading order 
in the hard scale, the SSA is given by 



A **&h-M = (_ 2 ) x 
£ Q e 2 /d 2 P /l± /d 2 p T /d 2 K T sin( 



(2) 



) sin(^ PT - fa) feL S^(zp T + K T - P hx ) xf#(x, p 2 T ) D<£(z, K 2 T ) 



J2 a e 2 a xf?(x)D<t(z) 

where we cancelled out the Q 2 and y-dependent factors (which describe the unpolarized partonic subprocess and are 
the same in the numerator and denominator), and Pr denotes the azimuthal angle around the z-axis e3 of the parton 
struck from the target nucleoli. In a full description of the process the transverse parton momenta in f^f(x, p 2 ") and 
D^(z, K 2 ^) would be convoluted with the soft factors [12-14] instead of the simplifying 5- function. 

However, this simplification is not yet sufficient for the purpose of extracting the Sivers function. In the numerator 
on the right-hand-side of Eq. (2) the integrals convoluting the transverse momenta cannot be solved, unless one knows 
fxr(x, p|i) and Df(z, K 2 ^), which is not the case. The situation would be different if in the SSA in Eq. (1) in addition 
had been introduced a power of the transverse hadron momentum |PhJ_|- 

The leading order expression for such a weighted SSA is given by [15] 
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where the transverse moment of the Sivers function is defined as 
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Preliminary HERMES data on the SSA weighted in this way were pre- 
sented in [8] . From the point of view of minimizing the model-dependence 
in the analysis [32] it is preferable to consider data on an appropriately 
transverse- momentum- weighted SSA [15], such as in Eq. (3). However, a 
wealth of data on the SSA in SIDIS have been analyzed and presented 
without the "appropriate transverse momentum weights" [4-7, 9-11]. 



FIG. 1: Kinematics of the SIDIS process 
Ip —* I'hX and the dehnitions of the azimuthal 
angies in the lab frame. 
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Many models for the transverse momentum dependence of distribution and fragmentation functions were considered 
in literature [26-31, 33-35]. Among the most popular models is the Gaussian ansatz, which has two important virtues. 
It describes successfully the distributions of low (with respect to the relevant hard scale Q) transverse momenta in 
various hard reactions — see, for example, Rcf. [23]. It also allows us to perform analytically the integrals over 
transverse momentum. In the Gaussian ansatz one assumes that the transverse momentum and x- or z-dependence 
of distribution or fragmentation functions factorize, and that the distributions of the transverse parton momenta are 
Gaussian: 



,2^ _ ta(„\ rT /funpJ 



exp(-K2/i^2 ) 

D?(z,K 2 T ) = Dt(z) \£l Dl ■ (5) 

The Gaussian widths p 2 np and are also referred to as mean square transverse momenta of the unpolarized 
distribution and fragmentation functions, respectively, since for example 

/ 2 \ _ Jd 2 pr PTf?(x,p 2 T ) Gauss 2 , , 

(PrJunp- J d2pr fa {Xjp 2 T) ~ Pun P - W 

In general, the mean square transverse momenta could be flavour and x- or z-dcpcndcnt, a possibility that we will 
disregard. For later convenience let us also introduce the notion of a mean parton transverse momentum generally 
defined, and in the Gaussian model given as follows: 

, > = Jd 2 P T |prj /l(s,Pr) Gauss \/7T 

^ T/unp - f d 2 p T f?(x,p 2 T ) _ 2 PuDp - [<) 
Under the above assumptions the expression for the Sivers SSA in Eqs. (1, 2) is given by [29] 

4 si„(0-0 s ) _ , 9 ^ acauss Eg eg xf^ 1)a {x) Df(z) M N 



Here we ignore the resolution cuts applied by the experiments: Pf t ± > 50 McV at HERMES [5, 9], and Ph± > 100 MeV 
at COMPASS [10]. Taking such cuts into account would yield the same expression (8), however, with an aGauss given 
in terms of incomplete T-functions. The error introduced here by neglecting these cuts, however, is marginal as we 
shall estimate below. In the Gaussian ansatz (5) the transverse moment of the Sivers function (4) is given by 



2M% 



The reason why in Eq. (8) we prefer to work with f^ V>a {x) instead of f^f(x) will become clear later. Before we start 
to extract the Sivers function from the HERMES data [9] , it is necessary to fix or constrain the free parameters pg iv 
and Kp , preferably by (other) HERMES data for sake of consistency. The next Section is devoted to this task. 



III. UNPOLARIZED SIDIS, POSITIVITY CONSTRAINTS, AND THE LARGE iV c -LIMIT 

Let us assume the distribution of transverse momenta in /f (x, p 2 -) and D"(z, K 2 -,) to be Gaussian according to (5), 
and let us furthermore assume the corresponding Gaussian widths to be flavour and x- or z-independent. Then the 
average transverse momentum of the produced hadrons as function of z is given by 

(Ph±(*)) = ^y^'pLv + k^. (10) 

Fig. 2 shows the HERMES data on (Ph±(z)} for h — pions from Rcf. [6]. Strictly speaking these data were taken from 
a deuterium target, but we will ignore this fact (i.e., we neglect nuclear binding effects and use isospin symmetry). 
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Also we will ignore the fact that these mean values are not corrected for acceptance effects. What is important for us 
is that these data allow to fix the free parameters p 2 and Kp . A best fit yields 



P 2 unp = 0.33 GcV 2 , 
= 0.16 GcV 2 , 



(11) 



and is shown in Fig. (2) by a dashed line. We observe a good and for our purposes sufficient agreement. It is important 
to stress that the agreement could be improved at the prize of introducing a z- and/or flavour dependent Gaussian 
width for the unpolarized fragmentation function, but we will refrain from doing so and stick to our simple picture. 

It is instructive to compare our result (11) to the values extracted in 
Rcf. [33] under certain assumptions from EMC data [42] on the so-called <p hl (z)> in GeV 



Cahn effect [43]. There, p 2 unp = 0.25 GeV 2 and K%. = 0.20 GeV 2 won- 



found. With these numbers one obtains a (Phj_(z)) which describes the 
HERMES data almost as well as the direct fit in Eq. (11) — Fig. (2). 
(Notice that in the formalism of Ref. [33] the expression (10) holds ap- 
proximately upon the neglect of terms of 0(k 2 _/Q 2 ).) Thus, we are lead 
to the encouraging conclusion that the Gaussian ansatz for ff(x, p§,) and 
D%(z, K|) is compatible with SIDIS data from HERMES [6] and that the 
ansatz and the numerical values for the Gaussian widths are supported 
qualitatively by the analysis [33] of EMC data [42] on the Cahn effect. 

Unlike in the case of the unpolarized fragmentation and distribution 
functions, it is not possible to fix the Gaussian width of the Sivers func- 
tion from SIDIS data. As the HERMES data [9] show a non-zero Sivers 
effect, of course, the parameter pg iv cannot be zero [72]. However, there is 
also a non-trivial upper bound for pg iv due to positivity conditions. The 
positivity bound for the Sivers distribution function reads [44] 
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If we demand the inequality (12) 
(5) at any value of x and for all 
sufficient condition must hold 



to be satisfied in the Gaussian model 
Pt|, then the following necessary and 
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2: The average transverse momentum 
(z)) of pions produced in SIDIS as mea- 
sured by HERMES from a deuterium target 
[6] vs. z. The dashed curve is the (Ph±(z)) in 
the Gaussian model (10) with the parameters 
as fixed here, see Eq. (11). The dotted curve 
is the (Ph±(z)) Gaussian model with the pa- 
rameters as obtained from a study of data on 
the Cahn effect [33]. 
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This means that the Gaussian width is bound from below and above as follows 
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In particular, if f^(x) ^ as the data tell us, then < R < 1 must be satisfied. Thus, we see that pg iv is restricted 
to the range 



< p§ iv < 0.33 GeV 2 . 



(15) 



It is worthwhile stressing that there is a bound for pg iv even if fyft (x) is very small. 

For later convenience let us derive from (12) bounds for the transverse moment of the Sivers function. From the 
Gaussian version (13) of the positivity bound we obtain immediately 
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Notice that from (12) one also can derive a model-independent bound as follows. Using 
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and the definition (7) we obtain 
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(17) 



In the derivation of this bound no use was made of any transverse momentum model. Therefore, it must be valid in 
any model. In fact, by evaluating (pr)unp in the Gaussian model (7), we find that the model bound (16) is stronger 
than (17), i.e. if the transverse moment of the Sivers functions satisfies (16) then it fulfills automatically also (17). 

When using the inequality (13) in our approach with the Gaussian widths assumed to be x- independent, it is 
understood that the ratio of the Sivers function (or its moment) to /f (x) is to be evaluated at that point in x, where 
it takes its maximal value. 

The present SIDIS data with their sizeable error bars [9, 10] do not constrain fits for f^f for the separate flavours 
a = u, d, u and d assuming the effects of heavier quarks to be negligible. In fact, in Ref. [33] where this has been 
attempted all fitted distributions but were found consistent with zero. In this situation it appears appealing to 
invoke additional theoretical constraints. In particular, here we will use predictions from the QCD limit of a large 
number of colours N c . 

In this limit the nucleon appears as N c quarks bound by a mean field [45], which exhibits certain spin-flavour 
symmetries [46]. By exploring these symmetry properties it was proven in a model independent way that in the 
large- N c limit [41] 

KAt" - fiT d )(x, P 2 T )\ » I (At" + ftf){x, p T )| , (18) 

S v ' S v ' 

= 0{Nl) = 0(N 2 C ) 

or, cquivalcntly, 

frr(x,p^) — — AtO^Pt) modulo l/N c corrections. (19) 

The relations (18, 19) are expected to be valid within their accuracy in the region of not too small and not too large 
x satisfying xN c — (D(N®). Similar relations hold for the Sivers antiquark distributions [73] . 
Inspired by the large- N c relation (19) we choose the following ansatz: 

xf^ 1)u (x) = -xf^ 1)d {x) = Ax b (1 - x) 5 , (20) 

and set f^ 1 ^{x), as well as the Sivers distributions of heavier quarks, to zero. The shape of the Sivers function at 

large x cannot be constrained by the data [9, 10]. The large- a; behaviour f^ V>q {x) oc (1 — x) b can be justified under 
certain assumptions — see [32] and references therein. However, one may consider it here as a mere model ansatz. 

Given the size of the error bars of the present data the ansatz (20) and the above assumptions are not too restrictive. 
This was exemplified in Ref. [32] in a study of preliminary HERMES data [8] on the transverse-momentum- weighted 
Sivers SSA (3). 



IV. EXTRACTION OF THE SIVERS FUNCTION FROM HERMES DATA 

In the previous Sections we fixed the parameter K% in Eq. (11), but all we have been able to do so far concerning 
Psiv was t° constrain this parameter to the vague range in Eq. (15). However, the parameter pg iv appears explicitly 
only in the "Gaussian factor" aQauss defined in Eq. (8). For illustrative purposes we evaluate this factor for the 
(z 2 ) (0.4) 2 of the HERMES experiment in the range (15) of the a priori possible value for pg iv , and we observe that 
it is strongly constrained to be in the range 

0.72 < a Gauss < 0.83 . (21) 

Thus, as long as one is interested in extracting from the HERMES data [9] the transverse moment f^ 1 ^ a {x) under 
the assumption of the Gaussian ansatz (5), then the result is only affected to at most ±10% by variations of the 
parameter j»g iv within the poorly constrained range given in (15). 

Note that as soon as we extract a result for f^ r>a {x), then from the inequality (16) we immediately obtain bounds 
for p^ iv which are stronger than those given in (15). We proceed as follows. 
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FIG. 3: a. The best fit and the respectively 1- and 2-o range for the parameters A and b in the ansatz (20) for the Sivers function. 
b. The dependence of the parameter A on the Gaussian width p| iv characterizing the transverse momentum distribution in the Sivers 
function in the Gaussian model (5). The parameter b is practically pg iv -independent. 




FIG. 4: a. The u-quark Sivers function xf^^ u (x) as function of x, as extracted from the HERMES data [9]. Shown are the best St, and 
its 1- and 2-cr regions, b. Here it is shown that the absolute value of the extracted Sivers function does not exceed half of the positivity 
bound in the Gaussian model in Eq. (16). 
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FIG. 5: The azimuthal SSA A s ^ h 0s) as function of x for charged pions as obtained from the fit (23) in comparison to the Rnal 
HERMES data [9]. 
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We use the parameterizations [50, 51] for /"(x) and D®(z) at a scale of 2.5 GeV 2 which corresponds to the average 
Q 2 in the HERMES experiment. We insert the large- N c motivated ansatz (20) in the expression for the x-dependcnt 
Sivers SSA in Eq. (8) 

A sin(4>-4>s) _ ,_ 9 x / dz QGaussQ) Eg e a xf^ 1)a {x) (z) 

where we consider the z-dependence of the Gauss factor acauss and integrate within the cuts [9] of the HERMES 
experiment, 0.2 < z < 0.7. We fix K£ according to Eq. (11) and choose for the parameter p| iv a value out of the 

range (15) of the a priori possible values. Then we check whether for the chosen p| iv the extracted /{^ . (x) does 
satisfy the inequality (16) within its 2-cr-uncertainty. We obtain the following result which refers to a scale of about 
2.5 GeV 2 : 

x.f^ 1)u (x) = -xf^ 1)d (x) = -(0.17... 0.18) x°- 66 (l - x) 5 , (23) 
with p| iv = (0.10... 0.32) GeV 2 . (24) 

Several comments are in order. The total x 2 is about 2.2, i.e., the x 2 per degree of freedom is about 0.3. 

In Fig. 3a the 1- and 2-cr regions of the parameters A and b are shown for pg iv = 0.20 GeV 2 , which is the central 
value of p| iv in the range (24). The dependence of the parameter b on pg iv is negligible. The response of the parameter 
A to variations of pg iv is nearly linear — Fig. 3b. 

The range (24) consists of those values of p| iv that are allowed by the general positivity inequality in the Gaussian 
model in Eq. (13). The data [9] do not allow us to constrain the parameter pg iv more accurately than that. However, 

F ±(i)a 

a nauiaiaLiuij uuaci vc nuw nine tuc lit icouil iui / 

Psw ^ S onr y poorly constrained in the range (24). 

In Fig. 4a we show the extracted u-quark Sivers function xf^j?~^ u {x). For the curve labeled as "best fit" we have 
chosen j)g iv = 0.20 GeV 2 which corresponds to the central value in the range (24) allowed by positivity requirements. 
For the minima (maxima) of the 1- and 2-cr regions we have chosen the minimal (maximal) value of pg iv in the range 
(24). Thus, the displayed error bands contain both, the statistical error of the HERMES data [9] and the uncertainty 
due to the poorly constrained Gaussian width pg iv of the Sivers function in Eq. (24). 

Notice that strictly speaking we neglected a low-P/jj_ cut in Eq. (8) and one may wonder how large is the error we 
introduced in this way. With the results we obtained (neglecting such cuts) in Eqs. (11, 24) we find that taking this 
cut into account in the HERMES kinematics would change our results by about 1% (and about 2% at COMPASS to 
be discussed below). Thus, in the Gaussian model the neglect of the \ow-Phj_ cut is justified. 

The absolute values of the extracted Sivers functions /{^"(x) = —xf^ V>d {x) are restricted by the upper (Gaussian 
model) bound (16) given numerically by 0.23 /"(x) or 0.23 f\(x). Since /i(x) is smaller than /"(x) the bound is 
stronger for the cZ-quark. In Fig. 4b we see that the extracted Sivers function well satisfies this bound, and we remark 
that it does not even exceed half of the general bound (17) within its 2-cr uncertainty. 

Finally, in Fig. 5 we compare the Sivers SSA obtained on the basis of our fit (23) to the HERMES data [9]. Of 
course, in the SSA the effects of the (small) uncertainty of the parameter A in (23) and the (sizeable) uncertainty of 
the parameter p| iv in (24) cancel. Notably, the 1-cr error band for the ir~ SSA is much narrower than for the ir + SSA. 
This means that the ir~ SSA is more sensitive to 1/N C corrections, i.e., to deviations from the ansatz (20). We will 
discuss this point in detail in the next Section. 

As an intermediate summary we conclude that the HERMES data [9] are well compatible with the large- N c predic- 
tions (19) for the Sivers function [41] and that the fit (23) satisfies the positivity bounds [44]. Remarkably, the sign 
of the extracted Sivers function in Eq. (23) is in agreement with the physical picture discussed in [52]. We remark, 
however, that model calculations of the Sivers function [53-55] show no tendency to exhibit the large- N c pattern (19). 



it is satisfactory to observe how little the fit result for ' a in Eq. (23) is affected by the fact that the parameter 



V. CROSS CHECKS: HERMES DATA ON z-DEPENDENCE & COMPASS DATA 

In our fitting procedure we did not use the HERMES data [9] on the z-dependence of the Sivers SSA. These data 

could have been used as an additional constraint for the integrals of x f^j, 1 ^ a (x) in the range 0.023 < x < 0.4, which 
corresponds to the cuts in the HERMES experiment. This would have helped to improve the significance of the fit, 
considering that only few x-data points are available. Instead, let us use these data here as a valuable cross check 
of our approach. As the z-shape of the SSA is dictated by the unpolarized fragmentation function D®(z) and the 
z-dcpendence of the Gaussian factor ctGauss in Eq. (8), this is not only a cross check for the extracted Sivers function 
(23), but it also tests the Gauss ansatz (5), the consistency of the choice of parameters (11, 24), and the large- N c 
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FIG. 6: A s [ jy Wh_0s) as function of z. The data are from the HERMES experiment [9]. The curves show the 1-cr variation of our large-N c 
constrained fit (23) for the Sivers function, and the impact of varying the Gaussian width p| iv of the Sivers function in the range (24). 
The z-dcpendent data were not used for the fit, i.e., the comparison serves as a cross check for our results. 



ansatz (19, 20) itself. In Fig. 6 we confront our fit result (23) with the z-dependent HERMES data on the Sivers SSA 
[9]. We observe a satisfactory agreement. Notice that the impact of the poorly constrained Gaussian width pg iv of 
the Sivers function (24) is marginal. 

The smallness of the tt" SSA can be explained as follows. Since fy^ = —frr m the large-7V c limit, the asymmetry 

A™ (<K ~ 0s) (tt~) oc (f£T f - 5^i av ), whcrc D[ m = Df^~ = Df- 71 ' = D^ w+ , etc. (D^ nf = = Df^ + , 

etc.), denotes the so-called favoured (unfavoured) fragmentation function. At any z we have the inequality 
D™ f (z) < D[ av (z). However, due to the weighting by the square of the quark electric charges, the effects of the 
smaller "unfavoured" and the larger "favoured" fragmentation function become comparable and tend to cancel each 
other. As a result the SSA for tt~ appears rather small, becoming zero around z = 0.56 (for the parameterization in 
[51] at Q 2 = 2.5 GeV 2 ). The tendency for cancellation persists with inclusion of the l/N c corrections, which, however, 
shift the position of the zero — see Sec. VI. 

Results from the COMPASS experiment [10] yield an equally important confirmation for the large iV c picture of 
the Sivers function. In the COMPASS experiment a solid polarized 6 LiD target [56, 57] was used. Neglecting, to 
a first approximation, nuclear binding effects and using isospin symmetry, we observe that for the deuterium Sivers 
distribution function we have f^jP ~ fyjP + f^ n ~ f" T + ff T , and analogously for d, q, etc. Thus, the deuterium 
target is sensitive to the flavour combination which is suppressed in the large- N c limit, see (18), and for which our 
ansatz (20) yields exactly zero. This is in agreement with the present COMPASS data which shows a Sivers effect 
from deuterium target compatible with zero within error bars [10]. 



VI. WHERE ARE SIVERS ANTIQUARKS & 1/N C CORRECTIONS? 

In the ansatz (20) we neglected the Sivers distributions for antiquarks and for the strange and heavier quarks. 
The neglect of strange and heavier quarks is probably a good assumption at the present stage of art. However, 
neglecting the Sivers antiquarks need not be such a good approximation. This can be seen from the unpolarized 
distribution functions, since in the region of x < 0.15 the distribution of unpolarized d-quarks reaches 25% and more 
of the unpolarized <i-quark distribution. In fact, this is precisely the x-region where HERMES [9] observes the most 
significant Sivers effect (for 7r + ) — see Fig. 5. 

In order to gain a rough idea of the possible uncertainty introduced by neglecting the Sivers antiquark distributions, 
we make two simple models for the Sivers u- and ^-distributions, while keeping the quark distributions at the value 
given by our fit result (23). Model I is that Sivers g-distributions are just ±25% of the corresponding Sivers quark 
distributions. This assumption may be an overestimate in the region of larger x. Therefore, in model II we set the 
ratio of each Sivers q- to Sivers g-distribution to be the same as for the unpolarized distributions. Thus, we will 
explore the effects of assuming each of the following models for the g-Sivers distributions: 

!0.25 = const model I , 
(f?+ff)(x) ( 25 ) 
Ul Jl ') model II . 
/r+/ w 
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A sta(Ms) (x) (a) A^^-^Vx) (b) 




0.1 0.2 0.3 0.4 x o 0.1 0.2 0.3 0.4 * 



FIG. 7: ^ h as a function of x. The data arc from the HERMES experiment [9]. The curves show the effect of modifying the 

result of large-Nc constrained fit (23) by using the Sivers antiquark distribution functions given in Eq. (25). The figures show that the 
data [9] have little sensitivity to the Sivers q-distributions, so that the neglect of the antiquark distributions in the fit was justified. 



Note that the particular ansatz for model II ensures compatibility with the large- N c limit, where (/" + fi)(x) <~ 
0(N%) 3> (/" — fi){x) <~ 0(N C ). Thus, our model Sivers antiquarks satisfy large-iV c relations analogous to (18, 19). 
We also automatically preserve the sum rule [64] (see also [32]), 

£ J Ax f^ 1)a {x) = . (26) 

a—g,u,d, ... 

In the large- N c limit the gluon Sivers distribution is suppressed with respect to the quark one [32]. Our models (25) 
satisfy the inequality corresponding to (17). Thus, being compatible with all theoretical constraints we are presently 
aware of, the models (25) are well suited for our purposes. 

In Fig. 7 we show the modest effect of the Sivers ^-distributions as a function of x. The effect is quite small for the 
7r + SSA, and well within the experimental errors. The effect is more pronounced for the 7r~ SSA, but still remains 
within the experimental errors. The effect of Sivers antiquarks on the ^-dependence of the Sivers SSA is less visible, 
and we refrain from showing analog plots. As here the x-dependence is integrated over, the entire effects amount 
of altering the overall normalization of the SSA without qualitative and with only small quantitative changes to the 
picture in Fig. 6. 

Thus, we conclude that even sizeable Sivers antiquark distributions, as modeled in Eq. (25), cannot be resolved 
within the error bars of the present data [9]. This justifies a posteriori the neglect of Sivers ^-distribution functions 
in our fit ansatz (20) here or in Ref. [32]. In this way we confirm also the observation made in Ref. [33]. There an 
attempt was made to extract /-j^? for the separate flavours a — u, d, u and d, and the Sivers-g distributions were 
found to be consistent with zero with large uncertainties. 

The difficulty to access antiquark distribution functions in SIDIS data is more general than the example encountered 
here. For example, in longitudinal double spin asymmetries All oc ^2 a e^g1(x)Di(z) it is presently [58] not possible 
to resolve a flavour asymmetry in the helicity sea as sizeable as predicted in the chiral quark-soliton model [59] , namely 
larger [60] than the known flavour asymmetry in the unpolarizcd sea — 9i)( x ) > |(/" — ff){x)\. As in the case of 
the much better known distribution functions f\(x) and (x) (see [61] for the specific example of helicity sea flavour 
asymmetry from the chiral quark-soliton model), it is necessary to study the Drell-Yan process in order to learn more 
about Sivers antiquark distribution functions. 

Next let us address the l/iV c -corrections. In our ansatz (20) we took literally the prediction from the large- N c limit 
[41] in Eq. (19), disregarding corrections which are generically of 0(1/N C ) w 30%. In order to have an idea of the 
effect of these corrections, let us assume that the flavour singlet Sivers distribution is not exactly zero but suppressed 
by exactly a factor of l/N c with respect to the flavour non-singlet combination according to Eq. (18). That is, 



(fir 1)u + f# 1)d )(*) = ± jfytr ,u - m L,a ){x) , m 

where we use f^^ u (x), /^^(x) from Eq. (23) and set N c = 3. Of course, it would be naive to expect that the 
different flavour combinations behave precisely as in Eq. (27). However, the scenario in Eq. (27) is by construction 
well suited to indicate the order of magnitude of the effect. In fact, when considered from such a qualitative point of 
view, large- N c relations for parton distributions are observed to be well satisfied in nature [62]. 



1 , r-L(l)u f -L(l)d\ 
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FIG. 8: a and b. A^ Wh 0s) as function of x. The data are from the HERMES experiment [9], The curves show our best Gt (23) for 
the Sivcrs function, and the effect of 1 / N c -corrections as modeled in Eq. (27). The figures show that the statistical uncertainty of the data 
[9] is of comparable magnitude as the 1/N C -corrections, c and d. The same as in Figs. 8a and 8b but for the z-dependence of the Sivers 
SSA. 
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FIG. 9: The Sivers SSA at COMPASS. The theoretical curves indicate the order of magnitude of the effect according to the rough model 
for 1/N C -corrections in Eq. (27). Thus, taking the largc-N c counting rule (18) into account we conclude that the HERMES [10] and 
COMPASS [10] data are compatible with each other. 



Note that (27) is compatible with positivity (17). However, in order to comply with the sum rule (26), we must 
introduce a gluon Sivers distribution function equal in magnitude but of opposite sign to the quark flavour singlet 
Sivers distribution. 

Fig. 8 shows the effect of l//V c -corrections as modeled in Eq. (27). The positive (negative) sign in Eq. (27) 
corresponds to the upper (lower) curves in Figs. 8a and 8b. What we obtained in this way is an "error band" of 
much the same size as the 1-cr uncertainty of our fit in Fig. 5. A look at the z-dependence of the Sivers SSA fully 
confirms these findings — see Figs. 8c and 8d. When the HERMES z-cuts are used, (aGauss-Di av ) ~ 2(aGauss-Di nf ) 
to a good accuracy with the parameterization [51]. If this were exactly the case, then for the minus sign in Eq. (27) 
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the 7T~ Sivers SSA would become exactly zero, as shown in Fig. 8b. More generally, we see that the vanishing of the 
7r~ asymmetry at some value of z is a feature which is robust against l/A^-corrections (and antiquark effects). The 
precise position of this zero, however, is very sensitive to l/7V c -corrections — Fig. 8d. 

On a deuteron target, the leading 1/N C prediction gives zero for the 
SSA, so that the 1/N C corrections are all that remain. Assuming for sim- 
plicity that the positive and negative hadrons identified at COMPASS are 
mainly pions, we obtain in our rough model (27) for l/A c -corrections the 
results shown in Fig. 9. The positive asymmetries are due to choosing the 
positive sign in Eq. (27). Clearly, we see that the COMPASS data [10] are 
compatible with the large- N c corrections being of a magnitude compatible 
with out model. 

Thus, the reason why our large-A c approach works here, is due to the 
fact that current precision of the first experimental data [9, 10] is com- 
parable to the theoretical accuracy of the large- N c relation (19). This is 
illustrated in a different way by Fig. 10 which shows that a variation of 
the best fit (23) of ±30% corresponding to the generic size of 1/N C cor- 
rections is of similar size and shape as the 1-a region of the fit. In other 
words, l/N c corrections (and antiquark effects) cannot be resolved within 
the error bars of the data [9, 10]. In future, with increasing precision of the 
data, it will certainly be necessary to refine the fit ansatz (20) to include 
1/N C corrections and antiquarks. In fact we have found that our fit (23) 
to the final HERMES data [9] is also compatible with the most recent and 
substantially more precise preliminary HERMES data given in [11]. 




(1 ± 30% jy. best fit 
1-G region of the best fit 
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The 1-a region of our best fit 



FIG. 

for ff^ 1 ^ u (x) in comparison to a variation of 
±30% of our best fit (23), which roughly corre- 
sponds to the generic size of 1 /N c corrections. 



VII. FURTHER TESTS OF THE GAUSSIAN ANSATZ IN SIDIS 



We have seen that the Gaussian model (5) for the distribution of the transverse parton momenta provides a 
satisfactory effective description of the HERMES data on the average transverse momentum of pions produced in 
SIDIS [6] . In order to increase our faith into the applicability of the Gaussian model and — what is equally important 
- to find out its limitations, it is necessary to make further tests using, e.g., the available HERMES data both on 
polarized and unpolarized SIDIS. Although the true prediction in QCD, or any other field theory, is that there is 
approximately a power-law fall off at large transverse momentum, it can well be that a Gaussian dependence can be 
a useful approximation for transverse momenta that are low with respect to the relevant hard scale. This, in fact, 
corresponds to the situation in the HERMES experiment [9] where (Phi) ~ 0.4 GcV < \/(Q 2 ) ~ 1.5 GeV. 

In this context it would be interesting to study the average transverse momentum square (P^±(z)) of the produced 
hadrons given by (actually this relation is of more general character [63] and manifestly valid in the Gauss ansatz) 



(pi ± ^)) = ^pL P + KE 1 . 

In particular, (Ph±(z)} and (P"h±( z )) are related to each other in the Gaussian model by 

(Ph±(z)) 2 G = ss \ (P£j_(*)) ■ 



(28) 



(29) 



To test the assumption of flavour independence, it would be useful to study (Ph±(z)) and (P^ ± (z)) separately for 7r+, 
7t°, 7r~ or the kaons. Finally, one could study the average transverse momentum of different hadrons averaged over 
z but as function of the respective x-bin, in order to test the assumption of an x- and flavour independent Gaussian 
width of the unpolarized distribution function. 

Given the unsatisfactory situation with the poorly constrained Gaussian width of the Sivers function, it is of 
importance to further constrain this parameter by means of data — in particular on the Sivers SSA as function of 
Ph±- Appropriate HERMES data for this as well as the following suggestion are in principle available [8, 9]. 

Let us define the following "Sivers-mean-transverse-momentum": 



(P h ±(z))sw = 



• h±ti sin(4 - <t> s ,i) [^v T (4; 4> s ,i) - A^(4; <t>s,i + *0] 



sin (<^ - M (</>*; <f>s,i) ~ Nl(<f>i; ^ + tt)] 



Gauss 2 



z2 P 2 si- 



+ K 2 Dl 



(30) 
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FIG. 11: The fits from Rcf. [32] of f^ 1)a (x) 
to the P h ±-wcightcd preliminary HERMES 
data [8] lie within the 1-cr region of the fit 
obtained here to the final (non-P^^-wcightcd) 
HERMES data [9]. 



One can read off the result in the Gaussian model by comparing Eqs. (3) 
and (8). Note that (Ph±(z))si v is of the same structure as its "unpolarized 
analog" (Ph±(z)) in Eq. (7) up to the coefficient 2/^/tt. From a simul- 
taneous analysis of the Sivers SSA weighted with and without a power 
of transverse hadron momentum (Ph±(z))siv could be determined with a 
relative accuracy comparable to that of the Sivers SSA. 

An indirect but important test of the Gaussian ansatz for the Sivers func- 
tion can be made on the basis of the analyses reported here and in [32]. In 
Ref. [32] the preliminary HERMES data [8] on the transverse-momentum- 
weighted Sivers SSA (3) were used in order to directly extract /^^"(ar) 
without resorting to any model for the transverse parton momenta. Here 
we have used the final HERMES data [9] to extract fi^ a (x) under the 
assumption of the Gaussian model (5). Thus, if the Gaussian model works, 
then the two extractions must yield the same result. 

In [32] the same set of assumptions (large-iV c for q, neglect of q, etc.) 
was used as here. However, instead of determining the l-a region, two 
different ansatze were explored: One ansatz is as in (20), and another is 
as in (20) but with fixed b = 1. In Fig. 11 we see that these fits from 
Rcf. [32] of f^ 1)a (x) to the P, l± -weightcd preliminary HERMES data [8] 
are compatible within 1-cr with the fit to the final HERMES data on the 
Sivers SSA weighted without a power of Phi. [9] . 

This observation indicates that the Gaussian ansatz for the Sivers function is compatible with the HERMES data 
within the statistical accuracy of the data [8, 9]. Given the preliminary status of the data [8] it is not possible to 
draw a stronger conclusion from this observation at the present stage. In fact, the HERMES Collaboration does not 
recommend [65] the use of the preliminary data [8] since they are not corrected for acceptance effects specific to the 
P/i_i_ -weighting, i.e., absent in the non-P/j^-weighted SSA [9]. However, our comparison in Fig. 11 indicates that the 
systematic error due to these effects is less dominant than the statistical uncertainty of the data [8, 9]. 

Another important test for the Gaussian model and for the fit result (23) could use the HERMES data (and possible 
future data from COMPASS and Jlab) on the ir° Sivers SSA. Preliminary data on the ir° SSA with large statistical 
uncertainties were shown in [8]. Taking into account the l-a uncertainty of the fit (23) we obtain the results shown 
in Fig. 12. Worthwhile commenting is the ^-dependence of the 7r°-asymmetry. Since the unpolarized fragmentation 
functions of light quarks and antiquarks into 7r° are the same (and since we neglect the effects of strange and heavier 

quarks), D"(z) completely cancel out from the SSA, so that ^ s \z,ir a ) oc aGauss(z). Thus the z-shape of the 

neutral pion SSA is entirely predicted by the Gaussian model. (This remains true even if one considers effects of Ph± 
cuts. The effect of the small resolution cut applied at HERMES or COMPASS [9, 10] is negligible in Fig. 12, see 
also Sec. IV.) A precise measurement of the z-dependence of the n° Sivers effect could therefore also help to test the 
Gaussian Ansatz. Alternatively, one may combine appropriately ir + and ir~ data and use isospin symmetry to arrive 
at the same information. This may yield results with better statistical accuracy in the HERMES experiment [69] . 
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FIG. 12: The Sivers SSA for neutral pions at HERMES as function of x and z, respectively, as predicted on the basis of the fit result 
(23). The error band arises from the l-a uncertainty of the fit (23). 
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VIII. CONCLUSIONS 

In this work we have extracted the transverse moment f^^ a of the Sivers function from SIDIS HERMES data [9] 
using a Gaussian model for the distribution of parton transverse momenta and employing predictions from the large- 
N c limit [41] as an additional constraint. We have shown that the Gaussian model provides a reasonable description 
of HERMES data on the transverse momentum distribution of the hadrons produced in unpolarized SIDIS. We 
constrained the free parameters of the Gaussian model consistently by HERMES data, which however, does leave the 
Gaussian width of the Sivers distribution poorly constrained. Nevertheless, the data [9] well constrain a fit of the 
transverse moment of the Sivers function. 

We have also shown that the HERMES and COMPASS data [9, 10] are compatible with each other and with 
predictions from the large- 7V C limit of QCD [41] within their statistical accuracy. We checked cxplictly that the effects 
of the simplifications we made are cither negligibly small, as for example the neglect of the experimental resolution 
cuts on the transverse hadron momenta, or well within the statistical accuracy of the data, as the usage of large- A^ c 
constraints or the neglect of Sivers antiquark distributions. We provided cross and consistency checks for the fit result 
by studying the z-dependence of the HERMES data on the Sivers SSA [9], and made suggestions how to further test 
the applicability of the Gaussian model in SIDIS. 

The main differences of our approach as compared to the similar works [33-35] are the use of the large- N c constraints, 
and the choice of a different model for transverse parton momenta and/or the way we fixed the respective parameters. 
Our fit is in qualitative agreement with extractions of the Sivers function [33—35] from the same [9] and from the more 
recent and more precise (but preliminary) HERMES data [11], see Ref. [70] for a detailed comparison. 

Of particular interest are studies of the Sivers effect in the Drell-Yan (DY) process, because the Sivers function (and 
other "time-reversal-odd" distributions) are expected to obey an unusual universality property, namely to appear with 
opposite signs in SIDIS and in DY [18] . The experimental check of this prediction is a crucial test for the understanding 
of the Sivers effect within QCD. On the basis of an analysis of the preliminary HERMES data [8] it was shown [32] 
that this change of sign of the Sivers function could be checked in p^ir~ — > l + l~X in the planned hadron-beam-mode 
of the COMPASS experiment [66], and in p^p — > l + l~ X in the proposed PAX experiment [67] (whose primary goal 
is to access the transversity distribution function in p^p^ — > l + l~X [68]). Our study of the final HERMES data [9] 
presented here confirms the analysis of Ref. [32] solidifying the conclusions made there — as do the estimates reported 
in Refs. [34, 35] obtained from analyses of the most recent preliminary HERMES data [11]. Estimates for COMPASS, 
PAX and RHIC on the basis of the results obtained here will be presented elsewhere [71]. 
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